Sputum quantity and adverse events during mechanical ventilation
Sputum and ventilator-associated complications
1 - Division of Infectious Diseases, Department of Medicine;
2 - Department of Biostatistics, Epidemiology, and Informatics;
Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA;
3 - Department of Medicine, University of Colorado, Aurora, CO;
4 - Division of Infectious Diseases, Mount Sinai School of Medicine, New York, NY
BJK is supported by the National Institute for Allergy and Infectious Diseases (K23 AI121485 and L30 AI120149), as well as Centers for Disease Control and Prevention (CDC) contract awards (BAA 200-2016-91964 and 200-2018-02919). BJK and DP are supported by the CDC Healthcare-Associated Infection Prevention Epicenters Program (U54CK000485).
The authors report no relevant disclosures.
Deidentified data, as well as analysis scripts and model code are available at github.com/bjkellyio.
ventilator-associated complications, ventilator-associated pneumonia, subglottic suction
Tracheal intubation and mechanical ventilation provide essential support for patients with respiratory failure, but the course of mechanical ventilation may be complicated by adverse events. The Center for Disease Control and Prevention (CDC)’s surveillance definitions of adverse ventilator-associated events (VAE) categorize events as infection-associated or not. We sought to understand how the quantity of sputum produced by mechanically ventilated patients relates to the onset of VAE and whether VAE are categorized as infection-associated or not.
We performed a case-crossover study including 87 patients with VAE and evaluated 848 days in the pre-VAE period at risk for VAE. Sputum quantity at time of VAE and in the immediate pre-VAE period was compared to sputum quantity in the preceding period. We used the number of as-needed subglottic suctioning events performed per calendar day to quantify sputum production. We used CDC surveillance definitions for VAE and to categorize whether events were infection-associated or not.
Sputum quantity is greater in the period immediately prior to VAE, compared to the preceding period. However, sputum quantity does not discriminate well between VAE that are categorized as infection-associated and VAE that are not.
Subglottic suction frequency may serve as a valuable marker of sputum quantity, which is associated with risk for VAE. However, our results require validation in a broader population of mechanically ventilated patients and intensive care settings.
Tracheal intubation and mechanical ventilation are life-saving interventions to address respiratory failure, but adverse events occur during the course of mechanical ventilation that can result in prolonged critical illness and increased risk of death. Traditional surveillance focused on infectious causes of adverse events, particularly ventilator-associated lower respiratory tract infection (VA-LRTI), including ventilator-associated pneumonia and tracheobronchitis. But other adverse events, such as congestive heart failure or septic shock, can also result in the need for increased ventilator support.
In 2013, the Centers for Disease Control and Prevention (CDC)’s National Healthcare Safety Network (NHSN) implemented a new surveillance paradigm for adverse ventilator-associated events (VAE), which advanced a broader, objective definition of ventilator-associated complications (VAC) (Klompas, 2013; Magill et al., 2013; Magill, Rhodes & Klompas, 2014). VAC are now defined as an increase of daily minimum positive end-expiratory pressure (PEEP) of at least 3cm H2O or an increase in daily minimum fraction of inspired oxygen (FiO2) by at least 20%, after at least two calendar days of stable or decreasing requirements. There are three tiers of VAC: (1) VAC without signs of infection; infection-associated VAC (IVAC), which are accompanied by leukocytosis/leukopenia or fever/hypothermia, as well as altered antibiotic treatment; and (3) possible ventilator-associated pneumonia (PVAP), which are IVAC with sputum studies suggestive of VA-LRTI. The NHSN definitions of VAE are well validated to predict ventilator-associated morbidity and mortality (Kobayashi et al., 2017; Meagher et al., 2019). VAE criteria also form the basis for promising novel antimicrobial stewardship strategies (Klompas et al., 2017).
However, current VAE surveillance definitions of VAE capture a different population than clinically defined VA-LRTI, which are typically diagnosed by physical examination and radiographic findings, as well as microbiologic test results (Klompas, 2014; Stevens et al., 2014; Fan et al., 2016; Ramirez-Estrada et al., 2018). More than 50% of traditionally-defined VA-LRTI fail to meet VAE criteria (Klompas, 2019). The best validated interventions to reduce VAE involve minimizing sedation and careful monitoring of volume status, rather than interventions specifically targeted to infection prevention (Lewis et al., 2014; Klompas et al., 2015; Klompas, 2015; Cocoros & Klompas, 2016; Vaewpanich et al., 2019; Liu et al., 2019; Rello et al., 2019).
Sputum production is a physical examination finding that is often used in the clinical adjudication of VA-LRTI, but its relationship to VAE defined by CDC NHSN criteria is unknown. It is challenging to quantify sputum production, but a marker of sputum production is the frequency with which subglottic suction is performed. Subglottic suction is a procedure performed on an as-needed basis by nurses, respiratory therapists, and physicians to clear lower airway secretions (Caroff et al., 2016). We sought to evaluate the relationship between the quantity of sputum produced and VAE. Specifically, we sought to define how sputum quantity at the time of VAE relates to sputum quantity before VAE, and to define whether sputum quantity discriminates between VAC and IVAC/VAP, using the frequency of subglottic suction events as a measure of sputum quantity. Here we report the results of a case-crossover study designed to address these questions.
We performed a case-crossover study to define the relationship between the quantity of sputum production (primary exposure), measured as the frequency with which nurses and respiratory therapists chose to perform as-needed subglottic suctioning, and the occurrence of VAE (primary outcome). We also analyzed whether VAE were categorized as VAC without suspected infection or as IVAC. For the primary anlysis, the day of VAE was the case period; the days prior to VAE were the control periods.
We evaluated 90 consecutive VAE that occurred across 87 unique subjects, over 3.5 years in the medical intensive care unit (MICU) of the Hospital of the University of Pennsylvania (HUP), a tertiary academic medical center. All patients who depend on mechanical ventilation at HUP were monitored for VAE, and all VAE were reported to NHSN. VAC and IVAC were captured by review of HUP NHSN reporting records. The study was approved by the University of Pennsylvania Institutional Review Board (protocol #824613).
Subjects were eligible for inclusion if they had been admitted to the HUP MICU, required mechanical ventilation for a period greater than two calendar days (i.e., at risk for VAE), and experienced a VAE event within the enrollment period.
VAE dates and VAC versus IVAC categorization were obtained from NSHN reporting records. Subject demographics, medical diagnoses, laboratory, and microbiology data were obtained from the Penn Data Store, a repository of electronic health record data. Subglottic suction events were extracted from daily flowsheets that report all activity related to mechanical ventilation.
The primary exposure was the quantity of sputum production, defined as number of as-needed subglottic suction events performed per calendar day. Calendar days were used because VAE (outcome of interest) are adjudicated per calendar day.
The primary outcome was VAE, and we compared the frequency of subglottic suction on the day of VAE and in the preceding 48 hours to the pre-VAE period. To evaluate whether quantity of sputum discriminates between VAC and IVAC, we similarly compared subglottic suction frequency across VAE categories.
We determined the necessary sample size by specifying the desired precision with which to detect the effect of sputum quantity to discriminate the pre-VAE and VAE risk periods (Rothman & Greenland, 2018). We anticipated that with 85 subjects we could detect an absolute risk difference attributable to sputum quantity of 15% with posterior certainty 4% - 26%. With the same subject number, we anticipated detecting a 25% absolute risk of IVAC versus VAC attributable to sputum quantity, with 95% posterior certainty 4 - 45%.
We organized data with R statistical software version 3.6.1 and generated plots using the “ggplot2” package (R Core Team, 2018). We evaluated potential confounders by comparing clinical variables of interest across exposure groups (high versus low frequency of subglottic suction at onset of mechanical ventilation) using Wilcoxon rank-sum testing (continuous variables) and Fisher’s exact test (categorical variables). For the primary analysis, we performed Bayesian multilevel logistic regression with random, subject-level slopes and intercepts. Like like conditional logistic regression, this approach accounts for subject-level strata in the case-crossover design (Maclure, 1991; Maclure & Mittleman, 2000; Rice, 2008; Hosmer, Lemeshow & Sturdivant, 2013). To permit partial pooling of parameter estimates and rigorous evaluation of model fit, Bayesian models were fit using Stan Hamiltonian Monte Carlo (HMC) version 2.19 via the “rethinking”, “rstan”, and “brms” packages (McElreath, 2016; Bürkner, 2017; Carpenter et al., 2017). Models were fit with 4 chains of 2000 iterations, then confirmed with HMC diagnostics (no divergent iterations, no iterations saturating maximum tree depth, Rhat statistic < 1.1 for all parameters, and E-BFMI > 0.2), and by examining the posterior distributions (Gabry et al., 2019). We characterized posterior credibility intervals and type S error (Gelman & Tuerlinckx, 2000).
Data, model code, and code used to produce the figures is available at github.com/bjkellyio.
We evaluated 90 consecutive VAE, which occurred across 87 unique subjects. For subjects with more than one VAE, we considered only the first VAE, so 87 VAE (47 VAC without associated infection; 40 IVAC, including 12 PVAP) were included in the final analysis. The median number of subglottic suction events at onset of mechanical ventilation (first full calendar day on ventilator) was 7. We first compared clinical characteristics across 39 subjects with high (greater than median) sputum quantity and 48 subjects with low (less than or equal to median) sputum quantity at the onset of mechanical ventilation (Table 1). We noted no significant differences in subject demographics, medical history, or features of critical illness at onset of mechanical ventilation. We did note a disparity in bacterial growth from respiratory culture in the 7-days prior to mechanical ventilation onset: only 1 (2.1%) of the subjects with low sputum quantity had recent bacterial growth from respiratory culture, versus 5 (13%) of subjects with high sputum quantity. But this difference was not statistically significant. We next evaluated the sputum quantity at onset of mechanical ventilation as a predictor of time to VAE. A Weibull parametric survival model showed no significant differences in time to VAE attributable to sputum quantity (Supplemental Figure 1).
Among the 87 included subjects, we evaluated 848 ventilator-days at risk for VAE (i.e., not including the first two days of mechanical ventilation). Though sputum quantity at onset of mechanical ventilation was not associated with time to VAE, we did observe an association between increased subglottic suction events and increased risk of both VAC and IVAC. Figure 1 shows the predicted probability of VAC or IVAC based on the number of observed subglottic suction events. With subglottic suctioning at 10 versus 5 events over a 24-hour period, the absolute probability of VAC increases by 3.07% (type S error < 0.024) the same day, 2.97% (type S error < 0.018) the next day, and 1.71% (95% CI includes 0) the day after. With subglottic suctioning at 10 versus 5 events over a 24-hour period, the absolute IVAC probability by 3.57% (type S error < 0.035) the same day, 3.43% (type S error < 0.018) increased risk of IVAC the next day, and 1.25% (95% CI includes 0) the day after. We also evaluated the impact of 15 versus 5 subglottic suction events over a 24-hour period and found VAC probability increases by 6.75% (type S error < 0.024) the same day, 6.67% (type S error < 0.018) the next day, and 3.67% (95% CI includes 0) two days later; with the same increase, IVAC probability increases by 8.1% (type S error < 0.035) the same day, 7.67% (95% CI includes 0) the next day, and 2.57% (95% CI includes 0) two days later. Given the small number of observed PVAP events, posterior credibility intervals were wide. But a positive relationship was observed between subglottic suction events and same-day PVAP risk, which increases by 2.82% (95% CI includes 0) as subglottic suction events increase from 5 to 10, and by 5.86% (95% CI includes 0) as events increase from 5 to 15.
Greater sputum quantity was associated with increased probability of VAC, IVAC, and PVAP. To determine whether sputum quantity could be used to discriminate IVAC from VAC, or PVAP from VAC/IVAC (not including PVAP), we evaluated a model comparing subglottic suction events in the period before IVAC versus VAC, and a model comparing subglottic suction events in the period before PVAP versus VAC/IVAC. No level of subglottic suction events reliably discriminated IVAC from VAC. Greater than 15 subglottic suction events was associated with greater probability of IVAC, but this threshold did not have high posterior certainty (Figure 2). Subglottic suction events did not discriminate PVAP from VAC/IVAC. In an exploratory secondary analysis, we evaluated a multivariable model including subglottic suctioning and other clinical variables that may contribute to respiratory failure to determine features not included in the IVAC definition (i.e., excluding temperature and white blood cell count) that discriminate IVAC from VAC. In our cohort of medical intensive care patients, we found that serum liver enzymes, alanine aminotransferase (ALT) and aspartate aminotransferase (AST), best discriminated IVAC and VAC, with elevated enzymes favoring VAC. The odds ratio of IVAC per standard deviation increase in ALT was 0.51 (95% CI 0.26 - 0.92), per standard deviation increase in AST was 0.56 (95% CI 0.26 - 0.996), and per standard deviation increase in the sum of AST and ALT was 0.54 (95% CI 0.27 - 0.98). Number of pressors and number of antibiotics favored VAC and IVAC, respectively. But unlike AST and ALT, these alone did not discriminate well (Supplemental Figure 2).
The observed association between sputum quantity and VAE risk prompted us to examine the relationship between sputum quantity and the behavior of clinicians. In a multivariable model utilizing all mechanical ventilation days (in contrast to the analyses above, which utilized only the pre-VAE period), we examined the relationship between subglottic suction events, fever, and the decision to order a sputum culture. Whereas fever was associated with the probability of a sputum culture order (1.9 (95% CI 1.5 - 2.5)), sputum quantity was not (Supplemental Figure 3).
Using a case-crossover design, we investigated the relationship between sputum production during mechanical ventilation and adverse ventilator-associated events (VAE) with the goal of understanding whether sputum quantity reflects underlying lung disease that will result in increased oxygen (FiO2 and PEEP) requirement. We quantified sputum production as the number of subglottic suction events per calendar day; we categorized VAE using the CDC’s NHSN definitions of VAC and IVAC. In this cohort of medical intensive care patients with VAE, sputum quantity was associated with both IVAC and VAC. Though the association with IVAC was stronger, sputum quantity did not discriminate well between infection-associated and non-infectious VAE. Secondary analyses revealed other clinical variables, prominently AST and ALT, that better discriminated IVAC from VAC in this cohort.
The relationship between sputum quantity and VAE is of interest because sputum quantity is often clinically interpreted as a sign of VA-LRTI, and the relationship between VAE and VA-LRTI remains unclear (Klompas, 2014; Stevens et al., 2014; Fan et al., 2016; Ramirez-Estrada et al., 2018). The finding that sputum quantity is higher at the time of VAE and immediately before VAE, relative to the earlier pre-VAE period, validates sputum quantity as a physical examination sign associated with risk for worsening respiratory failure. However, sputum quantity did not discriminate well between IVAC and VAC. Greater sputum quantity was associated with VAC without fever or leukocytosis, as well as with IVAC. There are several possible explanations for the discordance between sputum quantity and the NHSN criteria used to discriminate infection-associated adverse events. One possible explanation is that the number of subglottic suction events does not reliably represent the quantity of sputum produced; misclassification of sputum production may result in poor discrimination between IVAC and VAC. Another possible explanation is that NHSN IVAC criteria, which are known to exclude a large proportion of clinically diagnosed VA-LRTI (Klompas, 2019), exclude true VA-LRTI events that would drive a stronger association with sputum quantity. A third possible explanation is that sputum production, though commonly interpreted as a sign of infection, is nonspecific, or that sputum quality (e.g., purulence) is more specific than sputum quantity.
Several limitations of our study must be noted. The case-crossover design was chosen because VAE are rare outcomes, and the exposure of interest – sputum quantity – is dynamic and most likely to impact VAE risk acutely. However, our exclusive focus on VAE cases leaves open important questions about the generalizability of our findings across the broader population of patients who depend upon mechanical ventilation. Likewise, our cohort was enrolled from the medical intensive care unit of a single, academic medical center. Sputum quantity may differ significantly in surgical, neurological, or pediatric intensive care patients, and practices for management of tracheal secretions may differ across institutions. Our findings must be validated with studies of sputum quantity across a broader population of patients dependent on mechanical ventilation.
Nevertheless, our study demonstrates a significant, positive association between subglottic suction events and VAE risk, supporting the validity of subglottic suction frequency ascertained from electronic medical records as a means to quantify sputum production. The relationships between the respiratory tract microbiome, lung immune activation, and sputum production during mechanical ventilation warrant further study, as does the impact of sputum production on clinical behavior. The frequency of subglottic suction events appears to serve as a valuable marker of sputum quantity for these future analyses.
| Table 1: Sputum Production at Onset of Mechanical Ventilation | |||||
|---|---|---|---|---|---|
| Greater Than 7 Subglottic Suction Episodes in First Calendar Day | |||||
| Characteristic | Overall, N = 87 | N | FALSE, N = 481 | TRUE, N = 391 | p-value2 |
| age | 58 (48, 66) | 87 | 60 (52, 66) | 53 (44, 66) | 0.20 |
| gender | 87 | 0.57 | |||
| M | 45 (52%) | 23 (48%) | 22 (56%) | ||
| W | 42 (48%) | 25 (52%) | 17 (44%) | ||
| race | 87 | 0.35 | |||
| asian | 3 (3.4%) | 3 (6.2%) | 0 (0%) | ||
| black | 23 (26%) | 12 (25%) | 11 (28%) | ||
| unknown | 6 (6.9%) | 2 (4.2%) | 4 (10%) | ||
| white | 55 (63%) | 31 (65%) | 24 (62%) | ||
| copd | 1 (1.1%) | 87 | 0 (0%) | 1 (2.6%) | 0.45 |
| asthma | 13 (15%) | 87 | 10 (21%) | 3 (7.7%) | 0.16 |
| ild | 0 (0%) | 87 | 0 (0%) | 0 (0%) | |
| lymphleuk | 3 (3.4%) | 87 | 1 (2.1%) | 2 (5.1%) | 0.58 |
| dm | 11 (13%) | 87 | 4 (8.3%) | 7 (18%) | 0.21 |
| chf | 26 (30%) | 87 | 14 (29%) | 12 (31%) | >0.99 |
| cirrhosis | 12 (14%) | 87 | 8 (17%) | 4 (10%) | 0.58 |
| ss1 | 7.0 (4.0, 9.5) | 87 | 4.0 (2.0, 6.0) | 10.0 (8.0, 11.5) | <0.001 |
| vent_duration | 15 (9, 29) | 87 | 15 (9, 27) | 16 (10, 31) | 0.36 |
| vae_date_rel | 7 (4, 11) | 87 | 7 (4, 11) | 8 (5, 11) | 0.28 |
| ivac_yn | 40 (46%) | 87 | 19 (40%) | 21 (54%) | 0.27 |
| pvap_yn | 12 (14%) | 87 | 5 (10%) | 7 (18%) | 0.48 |
| fever_7before_vent | 15 (17%) | 87 | 9 (19%) | 6 (15%) | 0.90 |
| respcx_7before_vent | 9 (10%) | 87 | 4 (8.3%) | 5 (13%) | 0.51 |
| admit_wbc | 13 (9, 19) | 86 | 13 (9, 18) | 12 (9, 19) | 0.79 |
| Unknown | 1 | 0 | 1 | ||
| admit_cr | 1.22 (0.83, 2.08) | 87 | 1.13 (0.80, 2.01) | 1.33 (0.85, 2.14) | 0.57 |
| admit_ast | 41 (22, 74) | 59 | 42 (21, 110) | 41 (23, 53) | 0.44 |
| Unknown | 28 | 14 | 14 | ||
| admit_alt | 27 (16, 50) | 59 | 30 (18, 54) | 21 (15, 36) | 0.30 |
| Unknown | 28 | 14 | 14 | ||
| admit_fio2 | 100 (60, 100) | 87 | 100 (60, 100) | 100 (55, 100) | 0.76 |
| admit_peep | 7.50 (5.00, 10.00) | 86 | 7.25 (5.00, 10.00) | 7.50 (5.00, 10.00) | 0.67 |
| Unknown | 1 | 0 | 1 | ||
| norepi_after0_1 | 43 (49%) | 87 | 25 (52%) | 18 (46%) | 0.74 |
| epi_after0_1 | 9 (10%) | 87 | 5 (10%) | 4 (10%) | >0.99 |
| pit_after0_1 | 23 (26%) | 87 | 16 (33%) | 7 (18%) | 0.17 |
| dopa_after0_1 | 4 (4.6%) | 87 | 2 (4.2%) | 2 (5.1%) | >0.99 |
| vanco_before0_7 | 27 (31%) | 87 | 16 (33%) | 11 (28%) | 0.78 |
| vancopo_before0_7 | 1 (1.1%) | 87 | 0 (0%) | 1 (2.6%) | 0.45 |
| metro_before0_7 | 4 (4.6%) | 87 | 2 (4.2%) | 2 (5.1%) | >0.99 |
| cefaz_before0_7 | 1 (1.1%) | 87 | 0 (0%) | 1 (2.6%) | 0.45 |
| dapto_before0_7 | 0 (0%) | 87 | 0 (0%) | 0 (0%) | |
| linez_before0_7 | 2 (2.3%) | 87 | 2 (4.2%) | 0 (0%) | 0.50 |
| mero_before0_7 | 6 (6.9%) | 87 | 3 (6.2%) | 3 (7.7%) | >0.99 |
| piptaz_before0_7 | 13 (15%) | 87 | 7 (15%) | 6 (15%) | >0.99 |
| cefep_before0_7 | 12 (14%) | 87 | 5 (10%) | 7 (18%) | 0.48 |
|
1
Statistics presented: median (IQR); n (%)
2
Statistical tests performed: Wilcoxon rank-sum test; chi-square test of independence; Fisher's exact test
|
|||||
Panels depict the results of random effects binomial regression models relating VAE to the frequency of subglottic suctioning. Panels A, B, and C show the association between increased sputum quantity and absolute probability of VAC on the same day, VAC on the next day, and VAC two days later, respectively. Panels D, E, and F show the between IVAC and subglottic suctioning, and panels G, H, and I the association between PVAP and subglottic suctioning, in the same order. Best estimates of the association are represented with a dark line, surrounded by bands of posterior certainty.
The results of random effects binomial regression comparing the subglottic suctioning in the pre-IVAC period to the period before non-infectious VAC (panels A, B, and C), or the pre-PVAP period to the period before VAC and IVAC (excluding PVAP) (panels D, E, and F) are shown. Best estimates of the association are represented with a dark line, surrounded by bands of posterior certainty. Very frequent subglottic suctioning is associated with IVAC, but no number of observed suction events discriminates IVAC from VAC with high certainty.
Results of a parametric Weibull survival model, relating survival without VAE (days) to the number of subglottic suction events over the first full calendar day of mechanical ventilation, are shown, suggesting no association between admission sputum quantity and subsequent VAE. The best estimates of the association is represented with a dark line, surrounded by bands of posterior certainty.
Supplemental Figure 2. Exploratory analyses of clinical variables that discriminate IVAC and VAC suggest a role for AST, ALT, number of pressors, and number of antibiotics.
The association between the sum of AST and ALT and probability of IVAC versus VAC are shown, stratied by number of pressors and number of antibiotics. Best estimates of the associations are represented with a dark line, surrounded by bands of posterior certainty.
Results of a random effects binomial regression relating subglottic suction events and fever to the clinical decision to send sputum for bacterial culture are shown. Best estimates of the association are represented with a dark line, surrounded by bands of posterior certainty.
Bürkner P-C. 2017. Brms: An R package for bayesian multilevel models using stan. Journal of Statistical Software, Articles 80:1–28. DOI: 10.18637/jss.v080.i01.
Caroff DA., Li L., Muscedere J., Klompas M. 2016. Subglottic secretion drainage and objective outcomes: A systematic review and Meta-Analysis. Critical care medicine 44:830–840. DOI: 10.1097/CCM.0000000000001414.
Carpenter B., Gelman A., Hoffman MD., Lee D., Goodrich B., Betancourt M., Brubaker M., Guo J., Li P., Riddell A. 2017. Stan: A probabilistic programming language. Journal of statistical software 76:1–32. DOI: 10.18637/jss.v076.i01.
Cocoros NM., Klompas M. 2016. Ventilator-Associated events and their prevention. Infectious disease clinics of North America 30:887–908. DOI: 10.1016/j.idc.2016.07.002.
Fan Y., Gao F., Wu Y., Zhang J., Zhu M., Xiong L. 2016. Does ventilator-associated event surveillance detect ventilator-associated pneumonia in intensive care units? A systematic review and meta-analysis. Critical care / the Society of Critical Care Medicine 20:338. DOI: 10.1186/s13054-016-1506-z.
Gabry J., Simpson D., Vehtari A., Betancourt M., Gelman A. 2019. Visualization in bayesian workflow. Journal of the Royal Statistical Society: Series A (Statistics in Society) 182:389–402. DOI: 10.1111/rssa.12378.
Gelman A., Tuerlinckx F. 2000. Type S error rates for classical and bayesian single and multiple comparison procedures. Computational statistics 15:373–390. DOI: 10.1007/s001800000040.
Hosmer DW Jr., Lemeshow S., Sturdivant RX. 2013. Applied logistic regression. Wiley. DOI: 10.1002/9781118548387.ch9.
Klompas M. 2013. Complications of mechanical ventilation–the CDC’s new surveillance paradigm. The New England journal of medicine 368:1472–1475. DOI: 10.1056/NEJMp1300633.
Klompas M. 2014. Discordance between novel and traditional surveillance definitions for ventilator-associated pneumonia: Insights and opportunities to improve patient care. Infection control and hospital epidemiology: the official journal of the Society of Hospital Epidemiologists of America 35:1196–1198. DOI: 10.1086/677641.
Klompas M. 2015. Potential strategies to prevent ventilator-associated events. American journal of respiratory and critical care medicine 192:1420–1430. DOI: 10.1164/rccm.201506-1161CI.
Klompas M. 2019. Ventilator-Associated events: What they are and what they are not. Respiratory care 64:953–961. DOI: 10.4187/respcare.07059.
Klompas M., Anderson D., Trick W., Babcock H., Kerlin MP., Li L., Sinkowitz-Cochran R., Ely EW., Jernigan J., Magill S., Lyles R., O’Neil C., Kitch BT., Arrington E., Balas MC., Kleinman K., Bruce C., Lankiewicz J., Murphy MV., E Cox C., Lautenbach E., Sexton D., Fraser V., Weinstein RA., Platt R., CDC Prevention Epicenters. 2015. The preventability of ventilator-associated events. The CDC prevention epicenters wake up and breathe collaborative. American journal of respiratory and critical care medicine 191:292–301. DOI: 10.1164/rccm.201407-1394OC.
Klompas M., Li L., Menchaca JT., Gruber S., Centers for Disease Control and Prevention Epicenters Program. 2017. Ultra-Short-Course antibiotics for patients with suspected Ventilator-Associated pneumonia but minimal and stable ventilator settings. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America 64:870–876. DOI: 10.1093/cid/ciw870.
Kobayashi H., Uchino S., Takinami M., Uezono S. 2017. The impact of Ventilator-Associated events in critically ill subjects with prolonged mechanical ventilation. Respiratory care 62:1379–1386. DOI: 10.4187/respcare.05073.
Lewis SC., Li L., Murphy MV., Klompas M., CDC Prevention Epicenters. 2014. Risk factors for ventilator-associated events: A case-control multivariable analysis. Critical care medicine 42:1839–1848. DOI: 10.1097/CCM.0000000000000338.
Liu J., Zhang S., Chen J., Mao Y., Shao X., Li Y., Cao J., Zheng W., Zhang B., Zong Z. 2019. Risk factors for ventilator-associated events: A prospective cohort study. American journal of infection control 47:744–749. DOI: 10.1016/j.ajic.2018.09.032.
Maclure M. 1991. The case-crossover design: A method for studying transient effects on the risk of acute events. American journal of epidemiology 133:144–153. DOI: 10.1093/oxfordjournals.aje.a115853.
Maclure M., Mittleman MA. 2000. Should we use a case-crossover design? Annual review of public health 21:193–221. DOI: 10.1146/annurev.publhealth.21.1.193.
Magill SS., Klompas M., Balk R., Burns SM., Deutschman CS., Diekema D., Fridkin S., Greene L., Guh A., Gutterman D., Hammer B., Henderson D., Hess D., Hill NS., Horan T., Kollef M., Levy M., Septimus E., VanAntwerpen C., Wright D., Lipsett P. 2013. Developing a new, national approach to surveillance for ventilator-associated events*. Critical care medicine 41:2467–2475. DOI: 10.1097/CCM.0b013e3182a262db.
Magill SS., Rhodes B., Klompas M. 2014. Improving ventilator-associated event surveillance in the national healthcare safety network and addressing knowledge gaps: Update and review. Current opinion in infectious diseases 27:394–400. DOI: 10.1097/QCO.0000000000000083.
McElreath R. 2016. Statistical rethinking: A bayesian course with examples in R and stan. CRC Press/Taylor & Francis Group.
Meagher AD., Lind M., Senekjian L., Iwuchukwu C., Lynch JB., Cuschieri J., Robinson BRH. 2019. Ventilator-associated events, not ventilator-associated pneumonia, is associated with higher mortality in trauma patients. The journal of trauma and acute care surgery 87:307–314. DOI: 10.1097/TA.0000000000002294.
Ramirez-Estrada S., Peña-Lopez Y., Kalwaje Eshwara V., Rello J. 2018. Ventilator-associated events versus ventilator-associated respiratory infections-moving into a new paradigm or merging both concepts, instead? Annals of translational medicine 6:425. DOI: 10.21037/atm.2018.10.54.
R Core Team. 2018. R: A language and environment for statistical computing.
Rello J., Ramı'rez-Estrada S., Romero A., Arvaniti K., Koulenti D., Nseir S., Oztoprak N., Bouadma L., Vidaur L., Lagunes L., Peña-López Y., EUVAE Study Group. 2019. Factors associated with ventilator-associated events: An international multicenter prospective cohort study. European journal of clinical microbiology & infectious diseases: official publication of the European Society of Clinical Microbiology 38:1693–1699. DOI: 10.1007/s10096-019-03596-x.
Rice K. 2008. Equivalence between conditional and Random-Effects likelihoods for Pair-Matched Case-Control studies. Journal of the American Statistical Association 103:385–396.
Rothman KJ., Greenland S. 2018. Planning study size based on precision rather than power. Epidemiology 29:599–603. DOI: 10.1097/EDE.0000000000000876.
Stevens JP., Silva G., Gillis J., Novack V., Talmor D., Klompas M., Howell MD. 2014. Automated surveillance for ventilator-associated events. Chest 146:1612–1618. DOI: 10.1378/chest.13-2255.
Vaewpanich J., Akcan-Arikan A., Coss-Bu JA., Kennedy CE., Starke JR., Thammasitboon S. 2019. Fluid overload and kidney injury score as a predictor for Ventilator-Associated events. Frontiers in pediatrics 7:204. DOI: 10.3389/fped.2019.00204.
## R version 3.6.0 (2019-04-26)
## Platform: x86_64-apple-darwin15.6.0 (64-bit)
## Running under: macOS Mojave 10.14.6
##
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRlapack.dylib
##
## locale:
## [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] patchwork_1.0.0 brms_2.12.0 Rcpp_1.0.4.6 tidybayes_2.0.3
## [5] gtsummary_1.3.0 gt_0.2.0.5 forcats_0.5.0 stringr_1.4.0
## [9] dplyr_0.8.5 purrr_0.3.4 readr_1.3.1 tidyr_1.0.2
## [13] tibble_3.0.1 ggplot2_3.3.0 tidyverse_1.3.0
##
## loaded via a namespace (and not attached):
## [1] colorspace_1.4-1 ellipsis_0.3.0 ggridges_0.5.2
## [4] rsconnect_0.8.16 markdown_1.1 base64enc_0.1-3
## [7] fs_1.4.1 rstudioapi_0.11 rstan_2.19.3
## [10] svUnit_1.0.3 DT_0.13 fansi_0.4.1
## [13] mvtnorm_1.1-0 lubridate_1.7.8 xml2_1.3.1
## [16] bridgesampling_1.0-0 knitr_1.28 shinythemes_1.1.2
## [19] bayesplot_1.7.1 jsonlite_1.6.1 broom_0.5.6
## [22] dbplyr_1.4.3 shiny_1.4.0.2 compiler_3.6.0
## [25] httr_1.4.1 backports_1.1.6 assertthat_0.2.1
## [28] Matrix_1.2-18 fastmap_1.0.1 cli_2.0.2
## [31] later_1.0.0 prettyunits_1.1.1 htmltools_0.4.0
## [34] tools_3.6.0 igraph_1.2.5 coda_0.19-3
## [37] gtable_0.3.0 glue_1.4.0 reshape2_1.4.4
## [40] cellranger_1.1.0 vctrs_0.2.4 nlme_3.1-147
## [43] crosstalk_1.1.0.1 xfun_0.13 ps_1.3.2
## [46] rvest_0.3.5 mime_0.9 miniUI_0.1.1.1
## [49] lifecycle_0.2.0 gtools_3.8.2 zoo_1.8-7
## [52] scales_1.1.0 colourpicker_1.0 hms_0.5.3
## [55] promises_1.1.0 Brobdingnag_1.2-6 parallel_3.6.0
## [58] inline_0.3.15 shinystan_2.5.0 yaml_2.2.1
## [61] gridExtra_2.3 sass_0.2.0 StanHeaders_2.21.0-1
## [64] loo_2.2.0 stringi_1.4.6 dygraphs_1.1.1.6
## [67] checkmate_2.0.0 pkgbuild_1.0.6 commonmark_1.7
## [70] rlang_0.4.5 pkgconfig_2.0.3 matrixStats_0.56.0
## [73] evaluate_0.14 lattice_0.20-41 rstantools_2.0.0
## [76] htmlwidgets_1.5.1 processx_3.4.2 tidyselect_1.0.0
## [79] plyr_1.8.6 magrittr_1.5 R6_2.4.1
## [82] generics_0.0.2 DBI_1.1.0 pillar_1.4.3
## [85] haven_2.2.0 withr_2.2.0 xts_0.12-0
## [88] abind_1.4-5 modelr_0.1.6 crayon_1.3.4
## [91] arrayhelpers_1.1-0 utf8_1.1.4 rmarkdown_2.1
## [94] grid_3.6.0 readxl_1.3.1 callr_3.4.3
## [97] threejs_0.3.3 reprex_0.3.0 digest_0.6.25
## [100] xtable_1.8-4 httpuv_1.5.2 stats4_3.6.0
## [103] munsell_0.5.0 shinyjs_1.1